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TECHNICAL  REPORT  SUMMARY 

In  this  investigation,  we  are  regionalizing  the 
Arctic  region,  Siberia  and  the  Eurasian  continental  area 
using  seismic  surface  waves.  The  purpose  of  this  region- 
alization is  to  determine  the  structural  properties  of 
the  upper  few  hundred  kilometers  of  the  earth,  and  to  de- 
termine the  variation  of  these  properties  from  one  sub- 
region  to  another  in  the  area  under  investigation.  Once 
the  structural  properties  in  the  various  regions  have  boon 
obtained,  we  will  be  able  to  apply  our  (existing)  optimized 
computer  techniques  to  the  computation  of  accurate  theoret- 
ical seismograms  for  any  hypothetical  type  of  source  loca- 
ted anywhere  within  this  continental  area.  Those  theoret- 
ical seismograms  can  be  applied  directly  to  the  discrimin- 
ation problem  by  comparing  seismograms  computed  for  both 
earthquakes  and  underground  explosions  with  the  actual  re- 
corded seismogram. 

The  first  technical  problem  to  be  dealt  with, 
which  has  recently  been  completed,  is  the  computation  of  ex 
peri  mental  phase  velocity  curves  for  surface  waves  travers- 
ing all  of  the  regions  under  investigation.  We  are  now  en- 
gaged in  the  operation  on  the  complete  set  of  these  44  fun- 
damental Rayleigh  wave  dispersion  curves.  The  interpreta- 
tion of  these  data  will  give  us  the  regionalization  of  the 
area.  This  latter  step  is  now  being  performed  with  exist- 
ing computer  programs;  the  data  are  obtained  from  analysis 
of  seismograms  taken  by  WWSSN  instruments.  The  programs 


have  already  been  used  successfully  in  the  regionaliza- 
tion of  the  Pacific  Ocean  area  (Kausel,  1972;  Leeds, 

1973;  Kausel,  Leeds  and  Knopoff,  1974;  Leeds,  Knopoff 
and  Kausel,  1974). 

The  phase  velocity  curves  were  obtained  with  the 
single-station  phase  velocity  method,  which  is  described 
in  the  main  text  of  this  report.  Briefly,  the  method 
involves  scanning  our  microfilm  library  of  World  Wide 
Standardized  Sei smograph i c Network  (WWSSM)  records  for 
earthquakes  which:  (1)  occurred  within  , or  on  the  peri- 
meter of,  the  area  of  interest,  (2)  produced  good  long- 
period  surface  wave  records  at  WWSSN  stations  for  which  the 
epi center-to-stati on  lines  lie  within  the  regions  being  investi 
gated,  and  (3)  generated  good,  large  recordings  at  a suffi- 
cient number  of  stations  to  ensure  an  accurate  fault  plane 
solution.  When  a suitable  earthquake  has  been  found,  an  ex- 
tensive data  processing  and  data  reduction  system  is  applied 
to  transform  the  data  into  phase vel oci ty  curves  for  each  of 
the  selected  epicenter-to-station  lines  through  the  area  of 
i n teres t 

The  first  two  months  of  our  work  on  this  project  were 
spent  in  detailed  planning  of  our  attack  on  the  problem  and 
in  adapting  our  data  processing  and  reduction  system  to  the 
type  of  records  to  be  used  in  the  project.  The  main  points 
of  interest  here  are  given  at  the  end  of  the  final  section  of 
the  report. 

In  the  remainder  of  the  contract  period,  we  have  managed 
to  complete  all  of  the  data  analysis  (for  Rayleigh  waves 
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reccn  ded  by  the  WWSSN  instruments)  which  w ill  be  required  for 
the  regionalization  of  the  Arctic  region,  Siberia  and  the 
Eurasian  continental  area.  We  are  currently  involved  in  an 
intensive  program  of  inversion  of  these  data  *o  obtain  the 
desired  lateral  regionalization  of  these  areas.  In  the  final 
six  months  of  the  contract  period,  our  primary  goals  are: 
first,  to  complete  the  regionalization  of  Siberia  arid  the  Eu- 
rasian continental  area,  using  only  the  data  for  paths  limited 
to  these  regions;  second,  to  carry  out  the  regionalization  of 
the  Arctic  area,  using  only  the  data  for  paths  limited  to  this 
region;  and  third,  to  perform  a final  regionalization  of  Si- 
beria, the  Eurasian  continental  area  and  also  the  Arctic  re- 
gion, using  the  entire  set  of  data  we  have  available,  in- 
cluding many  profiles  which  include  both  continental  and 
oceanic  segments  between  epicenter  and  station.  Also,  by 
the  end  of  this  contract  period  , we  hope  tc  have  finished 
our  tests  of  the  records  from  the  new  high-gain  wide-hand 
instruments  at  CHG,  MAT,  EIL,  TOL  and  ICON.  If  these  tests 
are  satisfactory,  the  records  from  the  new  instruments  should 
provide  longer-period  information,  and  hence,  structural  infor- 
mation to  greater  depths. 

from  our  preliminary  regionalization  work  on  the  continental 
areas,  we  can  draw  a few  conclusions  already.  First,  we  note 
that  the  properties  of  much  of  the  Asiatic  upper  mantle,  av- 
eraged over  long  trans-Eurasian  distances,  arc  highly  consis- 
tent with  those  of  young  stable  regions  observed  elsewhere  (Knop- 
off,  1972;  Fouda,  1 973;  ICnopoff  and  I'ouda,  1974).  Second, 
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the  Tibetan  plateau  has  an  extremely  thick  oust,  perhaps  as 
great  as  100  km  or  more  from  surface  to  Moho.  Third,  the 
Alpine  folded  belt  of  Iran  and  Turkey  has  an  extraordinarily 
well -devel oped  low-vel oci ty  channel  in  the  mantle,  of  remark- 
able contrast  to  the  lid  above.  The  principal  problem  in  the 
inversion  has  been  the  construction  of  the  boundaries  to  the 
geological  provinces,  for  which  only  incomplete  information 
is  found  in  the  literature. 
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I.  Purpose  of  project 

We  are  studying  the  localized  properties  of  surface  wave 
propagation  in  the  Arctic  region,  Siberia  and  the  Eurasian  con- 
tinental area  using  regionalization  techniques  recently  devel- 
oped at  UCLA.  The  regionalization  yields  localized,  long- 

.-..period  surface-wave  dispersion  data  with  the  single- 

, station  method.  The  regions  of  high  seismicity  within 

and  around  the  region  to  be  investigated  are  quite  favorably 
located  for  the  application  of  the  single-station  technique. 

i 

The  dense  set  of  WWSSN  stations  located  around  the  perimeter 
of  the  region  ensures  that  we  can  acquire  sufficient  data 
for  the  application  of  our  regionalization  methods.  The 
new  high-gain  long-period  installations  in  Norway,  Spain, 
Israel,  Thailand,  and  Japan  arc  also  well  situated  for  the 
investigation.  The  long-period  data  from  these  latter- 
stations  may  provide  high  structural  resolution  at  depths 
. significantly  greater  than  were  previously  accessible  in 
surface  wave  studies.  Knowledge  of  the  regionalized  surface 
wave  and  structural  properties  will  allow  the  seismograph 
response,  to  an  arbitrary  source  within  the  region,  to  be 
synthesized  at  station  locations  on  the  perimeter. 


II.  Review  of  scientific  background 

Since  the  single-station  surface  wove  method  allows 
all  stations  to  he  located  at  the  edge  of  the  region  under 
investigation,  this  technique  is  ideal  for  the  proposed  study. 
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Brune  et  al  . (I960)  were  the  first  to  describe  the  single-sta- 
tion method.  Knopoff  and  Schwab  (1968)  corrected  and  extended 
this  description  to  take  into  account  the  frequency  depen- 
dence of  the  apparent  initial  phase  of  the  source.  Based  on  the 
Thomson  (1  950 ) -Has l;el  1 (1  953)  matrix  formulation  for  surface 

wave  propagation,  ilarkrider  (1  964)  developed  the  theory  for 
treating  the  surface  wave  response  to  buried  point-source 
singlets.  Ben-Menahein  and  ToksBz  ( 1 963)  then  developed  the 
formalism,  in  terms  of  the  singlet  response,  for  representing 
the  displacement  field  of  an  arbitary  force  system  in  a 
multi  lay ared  medium;  and  Ben-Henahem  and  Ilarkrider  (1  964) 
applied  this  formalism  to  obtain  the  surface  wave  response 
to  point-source  couples  and  double  couples.  Ilarkrider  (1970) 
later  made  certain  corrections  and  improvements  in  these 
last  results.  The  far-field  response  to  realistic,  dis- 
placement-dislocation faulting  was  shown  to  be  equivalent 
to  the  point-source,  double-couple  force  replacement  in  an 
unfaulted  medium  (Maruyama,  1 963;  Burridcje  and  Knopoff, 

1964),  and  the  single-station  method  was  on  firm  theoretical 
ground  for  flat,  multilayered  media.  By  means  of  tians- 
formation  techniques  (Biswas  and  Knopoff,  1970;  Schwab  and 
Knopoff,  1972),  it  is  possible  to  convert  point-source 
programs  for  a flat  structure,  into  those  useful  for  treat- 
ing sources  in  a spherical  structure  (Kausel  and  Schwab , 1 973 ) . 
These  programs  for  the  surface  wave  response  to  displacement 
dislocations  in  a radially  heterogeneous  sphere  are  now  a 
part  of  our  program  library.  We  are  using,  improving*  and 





optimizing  these  programs  as  part  of  the.  present  project. 

The  Initial,  large  scale,  lateral  regionalization 
of  a portion  of  the  earth  by  means  of  single-station  surface 

* 

wave  studies  was  performed  by  SantG  (1  960a;  1 9 C- 0 b ; 1961a; 

1961b)  in  ah  investigation  of  the  Pacific  Ocean  basin.  The 
bandwidth  of  his  recordings  and  his  data  processing  tech- 
niques limited  almost  all  of  his  dispersion  results  to  peri- 
ods less  than  about  40  seconds;  Santo's  regionalization  did 
not  yield  structural  information  below  the  crust  and  lid  of 
the  mantle.  Our  current  data  processing  techniques  (Biswas, 
1971),  and  the  availability  of  a nearly  complete  library  of 
long-period  seismograms  from  the  WWSSN  stations  have  made  it 
possible  to  extend  the  lateral  regionalization  to  a depth  o'f 
abqut  250  km.  The  first  application  of  our  single-station 
regionalization  techniques  was  an  investigation  of  the  upper 
mantle  structure  in  the  Pacific  Ocean  basin  (Leeds,  1973). 

The  results  of  the  work  of  Kausel  (1972)  showed 
an  apparently  continuous  gradation  of  Rayleigh-wave  phase  velo- 
cities from  low  values  observed  on  paths  close  to  the  Hast  Pa- 
cific Rise  to  relatively  high  values  on  paths  in  tiie  oldest  parts 
of  the  lithosphere.  Leeds  (1973)  increased  the  body  of  data  and 
performed  the  first  successful  inversion  of  the  data,  assuming 
that  regional  variations  in  structure  were  correlated  with  lith- 
ospheric age;  that  study  took  into  account  bathymetry  and  avail- 
able crustal  information.  The  number  of  degrees  of  freedom  in 
the  data  was  found  to  be  remarkably  small,  despite  the  large 


number  of  phase- vel oci ty  determinations.  The  inversion 
led  to  the  conclusion  that  the  data  sot  did  not  permit 
obtaining  detailed  information  concerning  the  bottom  of 
the  low-velocity  channel.  Furthermore,  the  lithosphere 
increased  in  thickness  monotoni cal ly  with  age:  at  the 
ridge  crest,  the  lithosphere  has  almost  zero  thickness; 
below  the  oldest  parts  of  the  oceanic  lithosphere,  the 
asthenosphero  has  almost  zero  thickness.  This  is  thus 
the  explanation  for  the  observation  that  the  phase  velo- 
cities change  systematically  with  distance  from  the  ridge 
crest.  The  model  of  lithospheric  thickness  as  a function 
of  spreading  age  is  consistent  with  a model  wherein  the 
1 i thosphere-asthenosphere  interface  is  at  the  solidus  for 
wet  peridotitc  if  the  thickness  of  the  slab  moving  lat- 
erally away  from  the  ridge  crest  is  relatively  thick,  of 

» 

a thickness  considerably  greater  than  the  depth  to  the  bot- 
tom of  the  asthenosphero.  Thus,  on  this  model,  there  is 
negligible  slip  between  the  asthenosphero  and  the  litho- 
sphere. Other  models  are  also  possible.  These  single- 
station  regionalization  techniques  arc  being  refined  for 
use  in  the  present  project. 

The  single-station,  surface  wave  regionalization  of 
an  area  is  based  on  the  assumption  of  great-circle 
propagation  paths,  and  the  assumption  that  the  phase 
travel  time  from  epicenter  to  station  is  just  the  sum- 
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of  the  travel  times  through  the  subdivisions  of  the  laterally 
heterogeneous  region  (Knopoff,  1969).  This  permits  the  region- 
alization to  be  expressed  at  each  frequency  as  a system  of  in- 
homogeneous equations;  for  each  path,  the  total  travel  time  is 
the  inhomogeneous  term,  the  distances  through  the  subdivisions 
are  the  coefficients,  and  the  slownesses  in  the  subdivisions  are 
the  unknowns.  Santfi  [equation'  (1),  1961b]  applied  this  technique 
to  the  Pacific,  using  group  travel  times  and  group  slownesses 
but,  due  to  considerations  involving  uniqueness  of  inversion  (Pi- 
la  n t and  Knopoff,  1970),  we  prefer  phase  travel  time  and  phase 
slownesses . The  solution  of  this  system  of  equations  yields  the 
experimental  phase  slownesses  associated  with  each  of  the  subdi- 
visions, in  principle;  however,  the  above  procedure  represents  a 
version  of  the  computer  program  which  is  not  pedagogical ly  accep- 
table. Instead,  we  prefer  to  consider  the  model  parameters  as 
the  primary  unknowns  in  the  inversion  and  to  use  the  phase  ve- 
locities therefrom. 

Techniq"es  for  obtaining  structural  information  from  disper- 
sion data  were  discussed  by  Knopoff  (1961,1962),  and  the  first 
programming  of  the  inversion  procedure  appears  to  have  been 
canied  out  by  Dorman  and  Ewing  ( 1 962  ).  In  these  papers, 
attempts  were  described  to  obtain  the  single  best  structure 
for  fitting  the  experimental  data.  Subsequent  efforts  were 
mainly  concerned  with  finding  the  set  of  structural  models 
which  fit  the  data  to  within  the  experimental  accuracy  (Keilis- 
Borok  and  Yanovskaya,  1967;  Press,  1968,1969).  This  work 


led  to  the  programming,  under  the  direction  of  Knopoff  and 
Keilis-Borok,  of  the  Hedgehog  inversion  package  (Biswas  and 

Knopoff,  1974;  Knopoff  and  Schluc,  1972).  This  tech- 
nique involves  seeking  structures,  which  satisfy  the  data 

to  within  the  experimental  accuracy,  by  means  of  a pointwise 
search  throughout  a multidimensional  paiameter  space.  This 
package  is  efficient  enough  for  routine  application  to  surface 
wave  dispersion  data.  Wo  have  now  introduced  variational 
parameters  into  the  procedure,  which  has  resulted  in  a 
significant  increase  in  computational  efficiency.  A descrip- 
tion of  the  ideas  involved  in  this  improvement  is  given  by 
Jackson  (1972);  the  technique  takes  into  account  the  redun- 
dancy in  the  data  as  well  as  numerical  instabilities  in  the 
solution.  • 

Upper  mantle  studies  in  Eurasia  appear  to  have  suffered 
due  to  the  Soviet  seismologists'  lack  of  efficient  long- 
period  instrumentation.  Also,  their  extensive  program  in 
deep  seismic  sounding  lias  undoubtedly  served  to  focus  interest 
upon  crustal  studies,  to  the  detriment  of  sub-Moho  investiga- 
tions. The  review  by  Arkhangel ' skaya  (1960)  discusses  the 
foreign  and  domestic  surface  wave  studies  performed  i.n 
Eurasia  prior  to  I960.  The  early  Soviet  studies  discussed 
in  this  review,  as  in  the  later  studies  to  be  mentioned  here, 
were  limited  to  short-period  invcstigations--usual ly  less  than 
40  seconds-- and  were  concerned  mainly  with  determining  crustal 
properties.  Of  the  early  Soviet  studies,  which  more  or -less 
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parol leled  the  work  in  the  west,  the  investigations  of  Sava  re ns ky 
and  Ragimov  (1  958;  1959),  Savarensky  ct  al_,  (1959),  and 
Savarensky  and  Sikharulidze  (1959)  clearly  demonstrate  this 
predominant  interest  in  crustal  features.  The  subsequent 
work  by  Popov  (1960),  Shechkov  (1961),  Savarensky  and 
Shechkov  (1961),  Shechkov  and  Solov'eva  (1961),  and  Shechkov 
(1964)  indicated  that  the  Soviet  surface  wave  investigations 
would  remain  focused  on  crustal  studies,  and  indeed,  this 
has  proven  to  be  the  case  up  through  the  most  recent  Soviet 
surface  wave  literature  available  to  us  (Savarensky  and 
Peshkov,  1968;  Sikharulidze  and  Makharadze,  1968;  Savarensky 
.et  al.,  1969;  Shechkov,  1970). 

All  of  these  crustal  studies,  as  well  as  Santo's  (1962; 

1965a;  1965b)  short-period  regionalization  efforts,.  will  be 
valuable  in  assisting  to  determine  which  shallow  structures 
to  combine  with  our  trial  parameteri zations  for  the  deep 
structures  in  the  inversion  portion  of  the  present  investiga- 
tion. 
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III.  Objectives,  methods  and  results  (of  the  first  18  months 
of  the  project) 

To  perform  the  regionalization  of  the  area  under  inves- 
tigation with  the  singl j-stati on  method,  we  must  know  the 
orientation  of  the  fault  planes  and  the  direction  of  dis- 
placement dislocation  at  the  focus.  These  provide  calcul- 
able corrections  due  to  apparent  initial  phase  of  the  Ray- 
leigh wave  signal.  Our  current  interpretation  procedures, 
for  obtaining  structural  parameters  as  a function  of  depth, 
are  based  solely  on  the  phase  velocity  dispersion  of  the  iso- 
lated fundamental  mode. 

Our  recent  large-scale  regionalization  work  (Biswas,  1971; 
Kausel  , 1 972  ; Leeds,  197V,  Fouda,  1 973  ; Kausel  et  al  . , 1974; 
Leeds  c_t_  aj_.  , 1974;  Knopoff  and  Fouda,  1 974  ) has  shown  the 
value  and  ease  of  handling  of  seismic  records  of  fundamental 
mode  Rayleigh  waves.  Our  initial  efforts  in  iho  present 
study  have  therefore  been  focused  on  this  well -tested  ap- 
proach . 

Only  the  long-period  records  from  the  47  Wl.'SSN  stations 
which  border  the  region  of  interest  have  been  used  thus  far 
in  the  study.  The  locations  of  these  stations  are  shown  in 
Fi gure  1 . 

In  order  to  obtain  the  desired  long-period  information, 
it  is  necessary  to  use  shocks  of  relatively  large  magnitude, 
yet  small  enough  to  be  on-scale  at  the  stations  of  the  net- 
work. ^n  example  of  the  large-magnitude  seismicity  of  the 
area  is  given  in  Figure  2.  The  epicenters  are  plotted  for 
earthquakes  having  magnitudes  between  5.9  and  G.G,  which 


occurred  during  the  interval  from  February,  1963  to  February, 
1967.  The  choice  of  this  range  of  magnitudes  is  governed  by 
two  considerations.  First,  experience  has  shown  that  good 
long-period  surface  wave  information  requires  events  with 
magnitudes  above  a certain  value  and,  of  course,  a shock 
which  is  so  large  as  to  send  the  instrument  off  scale  is  use- 
less for  our  purposes.  Second,  the  application  of  the  single 
station  method  requires  knowledge  of  the  focal  mechanism.  He 
must  therefore  use  events  large  enougli  to  allow  us  to  obtain 
an  accurate  fault  plane  solution  for  each  event  we  select  for 
processing.  In  addition  to  the  epicenters  shown  in  Figure  2, 
there  are  regions  of  high  seismicity  along  the  eastern  border 
of  the  Kamchatka  peninsula  and  along  the  Aleutian  arc,  which 
may  provide  useful  events  for  this  study. 

With  such  a dense  set  of  stations  around  the  area  to  be 
studied,  and  with  the  regions  of  high,  large-magnitude  seis- 
micity located  within  and  around  the  area,  there  has  been  no 
problem  in  obtaining  sufficient  data  for  the  project.  It  is 
interesting  to  note  that  the  area  is  almost  completely  en- 
circled: There  is  only  one  significant  gap  --  between  Japan 

and  Alaska  --  where  stations  do  not  exist.  However,  this  gap 
is  occupied  by  a region  of  high  seismicity  which,  in  the  sin- 
gle-station sense,  is  equivalent  to  having  a high  station 
density  in  the  region. 

There  are  far  too  many  cpicenter-to-stati on  paths  to 
show  individually,  but  the  limits  of  the  area  which  we  hope 
to  cover  finally  with  a dense  set  of  paths,  are  shown  in 
Figure  3 by  the  solid  lines.  The  solid  regions  are  those  of 


high,  large-magnitude  seismicity. 

The  first  event  processed  has  the  following  USC6S  speci- 
fication. 

March  7,  1966  - 21:29:17.4  GMT,  37.3°N,  114. 9°E, 
h = 3 3 km  , M=6 . 0 . 

Our  fault  plane  solution  for  this  event  is  given  in  Figure 
4.  This  first-motion  information  restricts  the  fault  plane 
solution  sufficient  well  so  that  only  this  information  is 
required  to  determine  the  required  focal  specification.  The 
location  of  the  epicenter  and  the  epi center-to-station  paths 
which  were  selected  for  processing,  are  given  in  Figure  5. 

For  all  of  these  recordings,  we  have  digitized  the  event  and 
impulse  response  from  the  seismogram,  isolated  the  fundamen- 
tal mode  with  frequency  filtering-time  windowing  techniques, 
corrected  for  the  instrument  response,  corrected  for  the  ap- 
parent initial  phase  due  to  the  source  mechanism,  and  compu- 
ted the  final  phase  velocity  dispersion  curves.  The  disper- 
sion data  for  these  eight  paths  are  in  final  form,  ready  to 
use  in  our  inversion  program. 

The  computation  of  the  correction  for  the  apparent  ini- 
tial phase  requires  the  specification  of  the  strike  of  the 
fault  plane  <p  , its  dip  6 and  the  direction  of  slip  A 
and  the  depth  to  the  focus  h . These  parameters  are  shown 
in  the  fault-plane  geometry  given  in  Figure  6.  For  the  above 


event,  our  fault-plane  solution  yields 
d)  = 122°  cast  of  north 
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and  inspection  of  the  Rayleigh-wave  amplitude  spectra  indi- 
cates a focal  depth  of 
h = 14.5  km 

The  second  event  selected  for  processing  has  the  follow- 
ing USCGS  specification. 

August  25,  1964  - 13:47:20.6  GMT,  78.2°N,  126. 6°E  , 
h=50  km,  M=6 . 1 . 

The  location  of  the  epicenter,  and  the  epicenter-to-station 
paths  which  were  selected  for  processing,  ore  given  in  Figure 
7.  The  fault-plane  solution  for  this  event  was  given  by 
Sykes  (1967)  as 

<i  = 4°  cast  of  north 
6=58°  west 

for  one  possible  fault  plane,  and 
<J>  = 22°  west  of  north 
6 = 54°east 

forthcothcr.  . 

Our  result  for  this  event  is 

<)  = 155‘'  east  of  north  • 

6 = 58° 

X * 260° 
h = 11.5  km . 

The  angles  <>  and  6 arc  well -cons trai ned  by.  ‘the  first- 
break  fault-plane  solution  given  in  Figure  8;  the  additional 

information  contained  in  the  Rayleigh-wave  amplitude  radio- 

% 

tion  patterns,  which  are  shown  -in  Figure  9,  is  required  to 
fix  X and  h . The  event  recordings  for  these  epicenter- 
to-station  paths,  and  the  corresponding  impul  se-respohs.o  re-- 
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cordings,  have  been  digitized;  the  fundamental  nude  has  been 
isolated  with  filtering-windowing  procedures,  instrument-re- 
sponse and  apparent-initial -phase  effects  have  been  removed, 
and  final  phase-velocity  dispersion  curves  have  been  prepared. 

. The  epicenter  and  propagation  paths  of  the  third  event 
selected: 

December  26,  1964  - 14:30:29.1  GMT,  51.8°N,  1^6. 8°E, 
h = 1 3 6 km,  M = 5.7  (USCGS) 

are  shown  in  Figure  10,  and  the  results  of  our  first-motion 
fault  plane  analysis  for  this  event  are  given  in  Figure  11. 

The  reduced  phase-velocity  curves  have  been  obtained  for  these 
paths  . 

The  epicenter  and  propagation  paths  for  the  fourth  event 
selected: 

March  31,  1969  - 07:15:54.4  GMT,  27.7°N,  34.0°E, 
h = 3 3 kin,  M=G  . 0 - 6.8 
are  given  in  Figure  12. 

Two  events  with  approximately  the  same  epicenter  were 
also  processed.  These  are 

February  6,  1965  - 01:40:33.2  GMT,  53.2°N,  161. 9°W, 
h = 3 3 km,  M = 6 . 4 - 6.7 

and 

February  6,  1965  - 16:50:29  GMT,  53.3°N,  161. 8°W, 
h = 33  km,  14-6 . 1 - 6.6 

The  paths  are  shown  in  Figure  13;  the  fault  plane  solutions 
are  given  in  Figure  14. 

Two  events  with  approximately  the  same  epicenter  were  the 
final  ones  processed.  The  USCGS  specifications 


are: 
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February  5,  1965  - 09:32:09.3  GMT,  52.3°N,  174. 3°F, 
h = 4 1 km,  M=5 . 9 - 6.5 
and 

February  6,  1965  - 04:02:53  GMT,  52.1°N.  175. 7°E, 
h=35  km,  M=5.9  - 6.0 

These  paths  are  shown  in  Figure  15;  the  fault  plane  solutions 
are  given  in  Figure  16. 

The  combined  Rayleigh  wave  propagation  paths  are  shown 
in  Figure  17,  where  the  continental  areas  are  emphasized,  and 
in  Figure  18,  where  the  Arctic  region  is  emphasized 
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Sample  phase  velocity  results  are  given  in  Figure 
19,  which  illustrate  the  variation  in  dispersion  for 
different  propagation  paths. 

Two  interesting  points  concerning  the  data- 
processing  techniques  have  come  to  light  during  the 
present  study.  Both  involve  the  accuracy  of  the  dig- 
itizations of  the  recorded  events.  First,  most  of  the 
event  records  we  have  used  arc  about  as  large  as  they 
could  be  without  going  off  scale.  (See  Figure  20  for 
an  example.)  This  has  necessitated  a change  in  otir 
data  processing  techniques  which  should  be  noted  for 
the  information  of  others  involved  in  this  type  of  work. 

In  the  past,  our  standard  procedure,  when  working  with 
small er-ampl i tude  recordings,  has  been  to  make  copies  of 
events  from  our  35mm  microfilm  records  of  the  WWSSN  seis- 
mograms. These  copies  are  made  with  a standard  microfilm 
reader-printer  (Itek  18.2/1  Reader-Printer),  and  the  events 
are  then  digitized  from  these  copies.  Wc  followed  this 
pi occduro  during  the  initial  phase  of  the  present  investi- 
gation, but  later  became  concerned  about  the  possibility 
of  distortion  in  the  copying  process.  Our  tests  comparing 
the  phase  velocity  results  obtained  from  full-size  record 
copies  provided  by  NOmA  with  the  results  obtained  from  our 
microfilm  copies  showed  this  concern  to  be  valid.  Our  con- 
clusion is  that,  when  working  with  large-amplitude  record- 
ings such  as  those  employed  in  the  present  study,  one  must 
use  full-size  record  copies;  one  must  not  use  a microfilm 
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copier  which  forces  the  analysis  to  use  several  prints, 
spliced  together  to  form  the  record  copy  from  which  the 
digitization  is  obtained. 

The  second  point  which  arose,  involving  accuracy 
of  digitization,  concerns  the  fact  that  the  direction  of 
swing  of  the  galvanometer  may  not  be  parallel  to  the  axis 
of  the  recording  drum  (James  and  Linde,  1971).  Although 
James  and  Linde  (1971)  term  this,  phenomenon  "a  source  of 
major  error  in  digital  analysis  of  WMSSN  seismograms,"  our 
tests  show  the  effect  to  be  negligible  on  the  phase  veloci- 
ties computed  using  the  single-station  method  for  epicenter- 
station  separations  of  a few  thousand  kilometers.  In  the 
case  of  the  poorest  galvanometer  alignment  we  encountered 
--  about  triple  the  normal  ramp  slope  of  0.3°  --  we  found 
only  negligible  differences  between  the  phase  velocity  curve 
obtained  with  the  correct,  digitization  base  line,  and  the 
curve  obtained  with  the  normal  ramp  as  a base  line. 
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FIG.  1 


FIG.  2 


FIG.  3. 


FiG.  4. 


FIG.  5. 


FIG.  G. 


FIGURE  CAPTIONS 

Locations  of  V.'WSSN  stations  and  new,  high-gain  installa- 
tions to  be  used  in  the  investigation. 


Example  of  the  large-magnitude  seismicity  of  the 
region  we  propose  to  study.  The  regions  of  high, 
large- magnitude  seismicity  along  the  eastern  border 

of  the  Kamchatka  Peninsula  and  along  the  Aleutian 

. 

arc  may  also  provide  useful  events  for  the  study. 

Limits  (solid  lines)  of  the  region  we  propose  to  cover 
with  a dense  set  of  epicenter-to-station  paths.  Solid 

regions  are  those  of  high,  large-magnitude  seismicity. 

• 

First-motion  information  and  fault  plane  solution 

for  event  occurring  at  21:29:17.4  GMT  on  March  7, 

' 

1 96G . 

' 

Location  of  epicenter  and  epicenter-to-station 
paths  selected  for  processing  from  event  occurring  at. 
21:29:17.4  GMT  on  March  7,  1906. 

Geometry  and  coordinate  systems  of  fault  plane, 
focus,  and  epicenter. 
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7.  Location  of  epicenter  and  cpi center- to-stati on 

p (\ t h s selected  for  processing  for  event  occurring 
at  13:47:20.6  GMT  on  August  25,  1964. 

8.  Ti rst-moti on  information  to  be  used  for  fault- 
plane  solution  for  event  occurring  at  13:47:20.6  GMT 
on  August  25 , 1 964 . 

9.  Rayleigh-wave  amplitude  information  to  be  used  for 
fault-plane  solution  inr  event  occurring  at  13:47:20.6 
GMT  on  August , 25 , 1 964  . The  central  set  of  radiation 
patterns  ore  the  results  of  theoretical  computations 
based  on  the  fault  plane  solution  given  in  the  text. 

The  other  four  radiation  patterns  depict  the  experi- 
mental results. 

10.  Location  of  epicenter  and  epicenter-to-station  paths 
selected  for  processing  from  event:  occurring  at 
14:30:29.1  GMT  on  December  26,  1964. 

11.  Results  from  first-motion  fault  plane  analysis  for 
event  occurring  at  14:30:29.1  GMT  on  December  26,  15C4. 

12.  Location  of  epicenter  and  epicenter-to-station  paths 
selected  for  processing  for  event  occurring  at 
07:15:54.4  GMT  on  March  31,  1969. 

13.  Location  of  epicenters  and  epi center- to-s ta* i on  paths 
selected  for  processing  for  events  occurrinj  at 
01:40:33.2  (dashed  lines)  and  16:50:29  GMT  (solid 
lines)  on  February  6,  1965. 
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FIG. 

FIG. 
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FIG. 


14.  Results  from  first-motion  fault  plane  analysis 
for  events  occurring  at  01:40:33.2  (Fig. 14a)  and 
16:50:29  GMT  (Fig. 14b)  on  February  6,  1965. 

15.  Location  of  epicenters  and  epicentcr-to-stati on 
paths  for  events  occurring  at  09:32:09.3  GMT  on 
February  5,  1965  (dashed  line)  and  at  04:02:53 
GMT  on  February  6,  1965  (solid  line). 

16.  Results  from  first-motion  fault  plane  analysis 
for  events  occurring  at  09:32:09.3  GMT  on  Febru- 
ary 5,  1965  and  at  04:02:53  GMT  on  February  6, 

1 965. 

17.  Combined  Rayleigh  wave  propagation  paths.  Con- 
tinental results  are  emphasized  here. 

18.  Combined  Rayleigh  wove  propagation  paths.  Arctic 
results  are  emphasized  here. 

19.  Sample  phase  velocity  results  for  the  paths  noted. 
The  paths  are  shown  in  Figures  5 and  7. 

20.  Example  of  large-event  recordings  used  in  the 
present  study.  This  is  the  recording  at  Kevo, 
Finland  (REV)  from  the  event  occurring  at  21:29:17.4 
GMT  on  March  7,  1966. 
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1 . 

Total  amount  of  present  award 

$114, 98G 

2 . 

Estimated  expenditures  and  commitments 

to  October  31 , 1974 
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Estimated  funds  required  to  complete 
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